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ABSTRACT 

Context. The generation and the amplification of magnetic fields in the current cosmological paradigm are still open questions. The standard 
theory is based on an early field generation by Biermann battery effects, possibly at the Epoch of Reionisation, followed by a long period of 
field amplification by galactic dynamos. The origin and the magnitude of the inter-galactic magnetic field is of primordial importance in this 
global picture, as it is considered to be the missing link between galactic magnetic fields and cluster magnetic fields on much larger scales. 
Aims. We are testing whether dwarf galaxies are good candidates to explain the enrichment of the Inter-Galactic Medium (IGM): after their 
discs form and trigger galactic dynamos, supernova feedback will launch strong winds, expelling magnetic field lines in the IGM. 
Methods. We have performed Magneto-Hydrodynamics simulations of an isolated dwarf galaxy, forming self-consistently inside a cooling 
halo. Using the RAMSES code, we have for the first time simulated the formation of a magnetised supernova-driven galactic outflow. This 
simulation is an important step towards a more realistic modelling using fully cosmological simulations. 

Results. Our simulations reproduce well the observed properties of magnetic fields in spiral galaxies. The formation and the evolution of our 
simulated disc lead to a strong magnetic field amplification: the magnetic field in the final wind bubble is one order of magnitude larger than 
the initial value. The magnetic field in the disc, essentially toroidal, is growing linearly with time as a consequence of differential rotation. 
Conclusions. We discuss the consequence of this simple mechanism on the cosmic evolution of the magnetic field: we propose a new scenario 
for the evolution of the magnetic field, with dwarf galaxies playing a key role in amplifying and ejecting magnetic energy in the IGM, resulting 
in what we call a "Cosmic Dynamo" that could contribute to the rather high field strengths observed in galaxies and clusters today. 

Key words, galaxies: formation - galaxies: evolution - galaxies: magnetic fields - methods: numerical 



1. Introduction 

The origin of magnetic fields in the Universe is a long-standing 
puzzle in cosmology. Although measuring magnetic field in 
cosmic structures is very challenging, we have strong evidence 
that magnetic fields exist in many galaxies, with magnitudes 
between 1 and 10 jiG at a l evel close to equipartition (see the 



review of Beck et al. Il996h . Magnetic fields of the order of 



sev eral tens of \iG are als o detected in large galax y clusters 
fsee lCarilli & Tavloill2002l or lGovoni & Ferettill2004l and refer- 
ences therein). However, in the very low density Intergalactic 
Medium (I GM), magnetic fi elds still remain undetected on 
large scale felasi et all 1999k . 

From a theoretical point of view, small-scale fluctuat- 
ing magnetic fields co uld be generated from qu antum effects 
during the Big-Bang (ITurner & Widrowl 1 19881 for example). 
The amplitude and the correlatio n length of such primordia l 
fields should both be quite small ( Grasso & RubinstelrfeOQll) . 
Magnetic fields greater than a few nG would also have been 



detected as a specific source of anisotropy in the Cosmic 
Microwave Background. More extreme scenarios can also be 
ruled out: for example, a uniform magnetic field greater than 
10" 1 jiG at our epoch would have led to unexpected anistropies 
in the expansion of the Universe (ICheng et al 11994b . Note also 
that the relic of any magnetic field generated in the early uni- 
verse, if it is greater than 10" 7 yuG today, would have broken 
the symmetry between left and right neutri nos (spin transition) 



during the pe riod of neutrino production (lEnqvist et al 
Sciamalll994l) . 



1993 
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When t he first structures col lapse (z - 10), the Biermann 
Battery (see lKulsrud et al, . 1997) is believed to have generated 
magnetic fields from pure collisional microscopic processes. 
These magnetic fields appear mainly at shocks and ionization 
fronts, where the motions of electrons and ions are decoupled, 
creating small microscopic currents. It was shown that this bat- 
tery effect can gen erate magnetic fields of the order of 10" 13 jjG 
field in the IGM dGnedin et al.ll2000h . More recently, shocks 
around primordial mini-halos and the associated first genera- 
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tion of massive stars were also co nsidered as imp ortant sources 
of magnetic field at high redshift dXu et al.ll20Q8h . 

In order to reconcile the very small values quoted above 
with observations, we need strong amplification processes up 
to the observed level of magnetic field strength in galaxies and 
clusters. Gravitational contraction of the frozen-in magnetic 
field will not be sufficient, since we expect in this case that 
B scales as p 2/3 . 

Using cosmological simulations of galaxy cluster for- 
mation, it was shown that one could reproduce the ob- 
served field magnitude by considering both gravitational 
contraction and turbulent stirring of the magnetised gas 
' Roettieer et lull 9991: bolag et al.ll 1 999L I2002L 120051: [siel et al.1 

et all 



2004 : lBruggenet al . 2005; S ubramanian et al.[l2006l:lAsai et al 



2007 : [Dubois & Teyssier! l2008al) . In cooling flow clusters, ra- 
diative losses are prob abily boosting this amplifi cation mecha- 



nism even further up dDubois & Tey ssier 
by recent observations dCarilli & Tayloi 



2008a) as confirmed 



2002). When an ac- 



tive galactic nuclei is releasing energy into the intra-cluster 
medium, a non-cool core develops and the magnetic field 
within the core is amplified by the jet-induced turbulence 
(Dub ois et al. I l2009l) . The conclusion of this series of paper is 
that one needs an initial magnetic field of at least 10" 5 yuG in 
the IGM to reach fiG values in cluster cores. Note that the fi- 
nite resolution of these numerical experiments leads probably 
to an underestimation of the turbulent amplification of the mag- 
netic field. Nevertheless, this justifies indirectly rather low val- 
ues for the magnetic field in the IGM, but significantly larger 
than Biermann Battery generated fields. 

In orde r to explain th e origin of this diffuse cosmic mag- 
netic field, iBertone et al. I d2006h proposed galactic winds as 
a possible solution to explain the enrichment of the IGM by 
metals and magnetic fields. In this scenario, magnetic fields 
are ampli fied inside gal actic discs by the so-called Galactic 
Dynamo dKulsrudlll999l) . and in a second phase, supernova- 
driven winds are expelling field lines into the surr ounding IGM 



The basic idea of the Galactic Dynamo ([Pa rker 



Ferrierd jl992al: iBrandenburg et alJ 1 19951: iKulsrudl 



1971 



1999; 



Shukurovll200l is that small-scale turbulence (precisely cy 



clonic motions that could result from sup ernova explosions, 
dsara et all |2004 iGissinger etaD l2QQ9h in the Interstellar 
Medium (ISM) modify sufficiently the small-scale magnetic 
field to make the large-scale component grow significantly. 
An additional term V x (aB) is added to the induction equa- 
tion that represents this dynamo amplification of the field, 
where a oc < u.V x u > tensor stands for the cyclonic 
motions of the gas . If one can correctly approximate the a 
parameter (IFerriere 11992A it is possible to show that the 
magnetic fi eld can grow in a few Gyr up to its equiparti- 
tion value (Ferriere & Schmitt l2000h . Note that the dynamo 
theory works only if the galaxy ca n loose magnetic helicity 
( Brandenburg & Subramanianll2005 ). even though it is still de- 
bated thatthe small scale saturation of the a effect really occurs 
( Fieldlll995l) . Numeric al simulations of a ve ry small patch (kpc 
size) of a galaxy from lGressel et alJ (120081) demonstrated that 
such a fast dynamo occurs if supernovae explosions can carry 
some magnetic helicity from the disc to the hot wind. However 
one has to reach a tremendous resolution power (parsec scale) 



in cosmological situations to provide such a long-term dynamo. 
Galactic winds have been considered as a nice explanation for 
suppressing some magnetic helicity from the disc, but no clear 
demonstration of the effect has been performed in a cosmo- 
logical context. Another issue with galactic dynamo is that 
the amplification mechanism is probably too slow, especially 
if one considers recent magn etic field observati ons in high- 
redshift galaxies beyond z = 1 dBernet et al.ll2008l) . These mea- 
surements suggest that magnetic amplification must occur very 
quickly, in less than 5 Gyr. It is also plausible that a fast growth 
of the field due to turbulent amplification in the first halos at 
redshift z > 10 cou ld provide an equipart ition field before the 
first galaxies form dArshakian et al.l l 20091) . Of course it has to 



be proven that small scale intense magnetic fields are able to 
generate this large-scale magnetic field observed in galaxies, 
for this purpose, a lternative theories must be explored. 



Parker! (119921) first suggested that the production of cos- 



mic rays in supernova remnants could generate stronger buoy- 
ancy in the ISM and amplify the magnetic field. Simulations 
of a cosmic ray pressurised medium have shown that it leads 
to a rapid growth of magnetic lines if they a re strongly dif- 
fused in the ISM dHanasz etal.ll2004 l2009allbh . iReesI d 1987k 
considered Biermann battery effects occuring at the surface of 
massive stars. More generally, one can consider any stellar dy- 
namo mec hanism as an effic ient amplification mechanism in- 
side stars dBrun et alJ 120041) . followed by the release of this 
magnetic field into the ISM, thanks to stellar winds or su- 
pernova explosions. Although the efficiency of stellar dy namo 
is rather uncertain dBrandenburg & Subramanian! 120051) . and 
the magneti c energy released into su pernova remants poorly 
constrained dKennel & Coronitill 19841 for the Crab nebula and 
iHelfand et al.ll200ll for the Vela nubula), this stellar origin rep- 
resents a very appealing perspective to explain a fast magnetic 
field generation in high-redshift galaxies. 

Galactic dynamos, or other field generation mechanisms in 
the ISM, together with galactic winds, should therefore play an 
important role in the generation of magnetic fields in the IGM. 
In this respect, this is probably inside dwarf galaxies that this 
double mechanism occurs: dwarf galaxies are very numerous 
in the early un iverse, and they host most of the cosmic material 
at early times dBertone et aDl2006l:lDonnert et al.ll2009h . Dwarf 
galaxies are also easily disrupted by galactic winds, in contrast 
to Milky Way like galaxies, fr om which metal and magn etic 
fields probably ne ver escaped ([ Dubois & Teyssier l2008bl and 
references therein). Kr onberg et al proposed a scenario 

where primeval galaxies launch strong starbursts that generate 
a 5. 10" 3 fiG IG magnetic field. They assume that galactic winds 
are able to reach Mpc scales and that they substantially amplify 
the field during the outburst. 

Understanding the evolution of magnetic fields in 
dwarf galaxies is therefore an important goal. Performing 
Magn etoHydrodynamics (M HD) numerical simulations with 
SPH, iKotarba et al] §009) studied the amplification of the 
magnetic field by differential rotation and by the spiral pat- 
tern of the gaseous disc. Another step was undertaken by 
Wang & Abell ( 20091) . who performed MHD numerical simu- 



lations with AMR of a dwarf galaxy formed by the collapse 
of a cooling halo. They studied in details the gas fragmenta- 
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tion and the magnetic field energy amplification due to turbu- 
lent, vortical motions induced by the gravitational instability in 
the disc. Starting with an initially uniform field of 10" 3 yuG in 
the halo, they reached equipartition after only a few rotations. 
They haven't simulated, however, the effect of su pernova feed- 



jpe 

back i n driving strong outflows out of the galaxy. iBertone et al 



d2006h . on the other hand, focused their study on the effect of 
galactic winds and how they might enrich the Universe with 
metals and magnetic fields. They used for that purpose a semi- 
analytical model, assuming equipartition magnetic fields inside 
their galactic discs and following the winds evolution through 
cosmic times. 

In this paper, our goal is to model the evolution of magnetic 
fields in a dwarf galaxy together with the formation of a galac- 
tic wind, in order to compute the amount of magnetic energy 
that escapes from the disc. This outgoing magnetic enegy flux 
is the key quantity that we need in order to estimate the effi- 
ciency of IGM enrichment. For that purpo se, we simulate an 
isolated , coo ling halo, in the same spirit of iDubois & Teyssier 
( 2008bl) and lWang& Abell d2009h . and form a small galactic 
disc in the halo centre. We solve numerically the ideal MHD 
equations, in the presence of self-gravity, with standard galaxy 
formation physics (cooling, heating, star formation, supernova 
feedback). The main parameter in our study is the starting mag- 
netic field that permeates the initial halo. We therefore vary its 
strength, assuming a simple but realistic initial field geometry. 
We will also consider the case of zero initial magnetic field in 
the halo, in a model where seed fields are introduced within su- 
pernova bubb les, impleme nting the stellar feed back scenario 
suggested by iReesI (1987) in a way similar to Hana sz et al 



(2009b|).InSect.[2]we present the numerical set-up of our simu- 
lations, as well as some details in our implementation of galaxy 
formation physics, and in Sect.[3j we discuss more specifically 
our initial conditions. In Sect. |4] and [5j we present our results 
on the magnetic amplification in galactic discs, and in Sect. 
we show how the IGM can be enriched by magnetised galac- 
tic winds. We finally comment the implications of this work in 
Sect.0 

2. Numerical methods 

2.1. Gas dynamics with magnetic fields 

We use the A daptive Mesh Refinement code RAMSES 
( Teyssierll2002l) to solve the full set of ideal MHD equations 



dp 
dt 



dpu 



+ V(pu®u-B®B + P tot I) 
dt 

dE 

— + V . ((£ + P tot )u - B(B.u)) 



dt 



dB_ 

~dt 



o, 


(1) 


0, 


(2) 


o, 


(3) 


o, 


(4) 



using a Godunov scheme with Constrained Tran sport that pre- 



serve s the divergence of the magnetic field (ITeyssier et al 
120061: iFromang et al.l l2Q06h . In the previous system, p is the 
plasma density, pu is the plasma momentum, B is the magnetic 



field, E = 0.5pu 2 + E± + B 2 /&7i is the total energy and E± is the 
thermal energy. This system of conservation laws is closed us- 
ing the Equation of State (EoS) for an ideal mono-atomic gas, 
where the total pressure is given by P tot = (y - l)E± + B 2 /8n, 
with y = 5/3. The MHD solver in the RAMSES code has been 
tested using ideal ized cases, as well as m ore realistic astrophys- 
ical situations in IFromang et al. I d2006h . It has been used for 
the first time in the cosmological context by Dubois & Teyssier 
( 2008al) to study the evolution of a cooling flow cluster, for 
which i t produced results in good agreement with previous 
studies (IDolag et al.l 120 05: Brii ggen et al. I l2005l) . Because it is 
based on Constrained Transport, the magnetic divergence ex- 
actly vanishes in an integral sense: the total magnetic flux 
across each AMR cell boundary is zero. The electric field is 
computed at cell edges, using a 2D Riemann solver, in or- 
der to upwind the edge-centered electric field with respect 
to the 4 neighboring cell states. The 2D Riemann solver is 
based on a generalisation of the ID HLLD Riemann solver of 
Miyoshi & Kusanol d2005l) . assuming a 5 -wave piecewise con- 
stant Riemann solution. We also account for self-gravity in the 
gas and stellar distribution, assuming a static potential for the 
dark matter component. The Poisson equation is solved with 
isolated boundary conditions, using Dirichlet boundary condi- 
tions given by a mutipole approximation of the mass distribu- 
tion in the box. The gravitational force is added as a source term 
in both the momentum and the energy conservation equations. 



2.2. Cooling, star formation and supernova feedback 

In order to describe the formation and evolution of a galac- 
tic disc, we model gas cooling in the halo as well as in the 
star forming disc using a look-up tabl e for the metallicity- 
depen dent radiative cooling function of ISutherland & Dopital 
(11993b for a mono-atomic gas of H and He with a standard 
metal mixture. This cooling function is added as a sink term in 
the energy equation ([3]). The metal mass fraction is advected 
as a passive scalar, with a source term coming from supernova 
explosions in the simulated ISM. 

Gas is not allowed to cool down to arbitraril y small tem- 
peratu re. Foll owing the general ideas prese nted in lYepes et al.1 
dl997l) and in lSpringel & Hernquistl d2003l) . we use a density- 
dependent temperature floor to account for the unresolved 
turbulence in the ISM, modelling this complex multiphase 
medium with an effective EoS. In this paper, we use the follow- 
ing temperature floor r m i n = 7o(ft/fto) 2/3 , where no is the star 
formation density threshold (see below) and To = 10 4 K. The 
effect of this temperature floor is to prevent the gas from frag- 
menting down to parsec scales, well below our r esolution limit. 
This i s in contrast with the study presented in IWang & Abell 
d2009h . where the gas temperature was allowed to cool down 
to 300 K, triggering strong disc fragmentation and massive 
clumps formation. 

Star formation is implemented using a standard approach 
based on the Schmidt law, spawning star pa rticles as a random 
Poisso n pr ocess. The reader can refer to iRasera & Teyssier 
d2006h and bubois & Tevssierl d2008bh for more details. We 
choose here the following star formation parameters: an effi- 
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ciency of 5% per free-fall-time and a star formation density 
threshold of no = 0.1 H. cm -3 . Each stellar particle represent 
a single stellar population of total mass m* = nomu^x 3 ^ 
10 4 M©. This stellar particle should therefore be considered as 
a Stellar Cluster (SC) rather than an individual star. 

After several Myr, the first massive stars explode in super- 
nova, giving rise to large expanding bubbles and ultimately to 
the formation of a large scale galactic wind. A mass fraction 
of ?7sn = 10% in massive stars M > 8 M Q (take n from any 



standard Initial Mass Function, e.g. lKroupall200ll) is consid 
ered for our stellar clusters. Modelling supernova feedback is a 
sensitive numerical issue that depe nds crucially on the adopted 
spatial resolution. It was shown by ICeverino & Klypinl (12009k 
that modelling feedback using a pure thermal dump of en- 
ergy works only at resolution smaller than 50 pc for a clumpy 
interstellar medium. In this paper, we used for our smallest 
cells a physical size of 150 pc. We therefore need to rely on 
a more sophisticated approach to model properly the forma- 
tion of supernova-driven bubbles. We adopt in this paper the 
now rather sta ndard kinetic feedback scheme, proposed by var- 
ious authors dS pringel & Hernq uistl 120031 Dubois & Teyssier 
2008bl: balla Vecchia & Schavell2008l) . For each stellar parti- 



cle formed, we also form in the parent gaseous cell an ad- 
ditional collision-less particle, mimicking the formation of 
a companion Molecular Cloud (MC). The mass of the MC 
is parametrized by m M c = W wtn*- The free param e ter tj^ , 
similar to the one defin ed i n ISpringel & Hernquistl ( 2003 ). 
iRasera & Tevssierl d2006h and balla Vecchia & Schavel J2OO8D . 
was choosen to be z/ w = 2. 

We assume here that after a delay of tsN = 10 Myr, super- 
nova feedback destroys the MC, and the corresponding mass, 
momentum, energy and metal content are released into the sur- 
rounding gas cells. Following Dubois & Tevssierl d2008bl) . we 
add to each fluid element a Sedov blast wave profile for den- 
sity, momentum and total energy, with a radius equal to a fixed 
physical scale of tsn = 300 pc, or 2 AMR cells. For metal en- 
richment, we assume a constant yield of vsn = 10% for any 
exploding star. A stellar cluster with mass m* = 1O 4 M and 
zero initial metallicity will release Vsn?7sn^* = 10 2 M© of met- 
als into the ISM. 

When we form a stellar cluster, we do not remove any 
large scale magnetic field component. We assume that at small 
scales, magnetic field is expelled from the collapsing MC by 
ambipolar diffusion, so that we can neglect the magnetic field 
that is trapped into the newly formed stars. Similarly, when 
supernova-driven blast waves form, we do not modify the mag- 
netic field in their surroundings. The magnetic field evolution 
is therefore driven only by the large scale velocity field. In only 
one case do we consider supernova seed fields, as explained 
in the Appendix, in o rder to test t he stellar magnetic feedback 
scenario proposed by lReesI (11987b . 



3. Initial and boundary conditions 

In this section, we decribe our initial conditions, especially our 
initial magnetic topology. We also decribe in details how we 
enforce the magnetic divergence free condition at the boundary 
of the computational domain. 



3.1. Initial hydrostatic halo structure 

We use the same set of initial co nditions as the one pre- 
sented in iDubois & Tevssierl ( 2008bl) . namely an isolated dark 
matter halo hosting a cooling gaseous component. This ide- 
alized situation have been used many times to study vari 



ous aspects of galaxy formation (Springel & Hern quistl I2003 



Dubois & Tevssierl 120085 iKaufmann et al.1 l2009h . Let us re- 
call here only the most important aspects. Dark matter is 
modelled with a sta tic density distribution following a NFW 
dNavarro et al.lll996h profile such that 



r/r s (l + r/r s ) 2 



(5) 



where p s and r s are the characteristic density and radius. This 
profile can also be parametrized in terms of concentration c = 
?*200 /r? and Virial mass as 



M 2 oo = 4np s r* 00 log (1 + c) + 



1 



1 +c 



(6) 



within the virial radius r2oo- The virial mass is defined here at 
an overdensity equal to 200 times the critical density p c . We 
also assume Ho = 70km/s/Mpc. 

The gas is initially in hydrostatic equilibrium with the same 
profile than dark matter, with a total gas fraction of ft = 15%. 
We consider that the halo is slowly rotating, with an angular 
momentum distribution consistent with the average specific an- 
gular momentum pro file measured in cosmological simulations 
dBullock etai]l200ll) 



j(r) = jn 



Af(< r) 
M 2 oo 



with a spin parameter 



J\E\ 



1/2 



GM 



5/2 
200 



(7) 



(8) 



Springel & Hernquistl d2003h 



Using this particula r settin 
and IDubois & Tevssierl (l2QQ8bh showed that large-scale out- 
flows are produced only in low-mass halos with M200 ~ 
10 11 M©. As we are interested in the role of galactic winds in 
the enrichment of the IGM, we consider in this paper only the 
case of a small mass halo with M200 = 10 10 M©. The halo spin 
parameter is chosen to be A = 0.1. 

The box size is set equal to 3 times the Virial radius or 
L = 150.7 kpc, and is covered by a coarse grid with 64 3 
cells. We consider 4 additional levels of refinement, so that 
our effective resolution reaches 1024 3 or Ax = 147 pc. Our 
refinement strategy is based on two criteria. First, we refine a 
cell if its mass exceeds 8 times our mass resolution indicator 
m res = 2 x 10" 6 M2oo- Second, we refine a cell is the local Jeans 
length falls below 4 cel ls, so that we prevent the disc from 



lengtn rails below 4 cel ls, so mat we prevent me disc trom 
artificially fragmenting dTruelove et aD 1 1 997b iBate & Burkert 
1997b lMachaceketal.ll200lk . The initial AMR grid was built 
from the coarse mesh with this refinement strategy. We start ini- 
tially with a total ~ 2.10 6 cells which number remains roughly 
constant with time. The number of cells in the maximum level 
of refinement (mainly located in the gaseous disc) is ~ 4.10 4 . 
After 3 Gyr the total number of stars is ~ 10 5 . 
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Fig. 1. (Oyz) slice of the initial magnetic field amplitude in log yuG units and magnetic vectors for #igm 
(left pannel) and t = 3 Gyr (right pannel). Picture size is 40 kpc. 



lO-^G at t = Gyr 



3.2. Initial magnetic field 

Although the choice of the initial magnetic field is quite crucial, 
the magnetic field topology is usually chosen quite arbitrarily 
in the literatur e. Using an init i al set up very similar to the one 
presented here JWang & Abelld2009h considered a constant ini- 
tial magnetic field B z - 10" 9 G. If the initial halo field is the 
product of the collapse of some large scale cosmological struc- 
ture, the field should rather scale as p^s • This scaling has been 
obser ved in cluster-size halo in cosmological MHD simula- 
tions dSigl et al.ll2004l; IBolag et alJEool iBriiggen et alll2005l: 
Dubois & Teyssierll2008al) . In the spherically symmetric set-up 
we are considering, each collapsing gas shell will compress 
and amplify significantly the frozen-in field, before any disc- 
like dynamo comes into play. The initial magnetic field am- 
plitu de profile is therefor e quite important. In a different con- 
text, ^ta^aetal] d2009l) considered the case of a pre-existing 
galactic disc. The initial magnetic field was also constant, but 
this time B x = constant, the field lines lying exacty in the disc 
plane. As argued by these authors, differential rotation quickly 
takes over and create a strong toroidal field in the star-forming 
disc, so that the initial topology is quickly forgotten. 

We nevertheless propose here to consider a more realistic 
initial magnetic field profile, with a dipole-like topology. In 
order to preserve the V.l? = constraint we compute the mag- 
netic field from its potential vector A 



B=VxA. 



(9) 



A is computed assuming that each cell in the computational 
domain contain a small magnetic dipole with strength propor- 
tional to Pg(s • Every dipole have the same direction along the z 



axis, so that they all add up to a global large scale dipole-like 
topology on the halo scale. The potential vector is written as 



A(x,y,z) = B ] 



IGM 



PgasM 


2/3 


(-y) 




+x 


. PlGM 




{ J 



(10) 



where pigm = ^bPc, with the baryon density Qb = 0-04 and r 
is the spherical radius. Each component of the potential vector 
is averaged over the corresponding cell edge, and the magnetic 
field is computed using the integral form of the curl operator on 
each cell face. One can see on Figure[T]that the initial magnetic 
field topology is roughly dipolar, preferentially aligned with the 
z axis, and that the field strength decreases strongly with radius. 

Figure [2] shows the magnetic field profile along the z axis 
for different values of the magnetic field in the IGM #igm- The 
magnetic field strength decreases down to 2 orders of magni- 
tude from the halo centre to the outer parts at 2r2oo- Results 
from cluster scale cosmological simulations seem to favor val- 
ues of B\qm ^ 10" 5 - 10" 4 //G in order to explain the observed 
magnetic field strength inside the core of large X-ray clusters 
dCarilli & Tavlorll2002l iGovoni & Ferettil 120041 and references 
therein). This value is however to be considered as an upper 
limit since these cosmological simulation have limited resolu- 
tion and are believed to underestimate the magnetic field am- 
plification from the initial IGM value to the final cluster's core 
strength. In this paper, #igm is a f ree parameter that we allow to 
vary between 10" 7 yuG, deep into the pure induction regime for 
which the Loren tz force remains negl igible and 10" 4 yuG, our 
maximum value. I Wang & Abelld2009l) initial magnetic field of 
B z = 10" 3 yuG in the halo corresponds roughly to our case with 
#igm - 3 x 10" 5 yuG, using the p 2/3 scaling and a mean over- 
density of 200 in the halo. 
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3.3. Zero-gradient boundary conditions 

A complex numerical issue for MHD flows is to design proper 
boundary conditions. We are dealing with an isolated halo, for 
which we consider that the environment is a very low density, 
homogeneous medium. It is convenient in this case to use out- 
flow boundary conditions: gas is allowed to flow freely out of 
the computational domain. Outflow boundary conditions can 
be implemented using a zero-gradient condition at the domain 
boundary for all hydro variable (density, pressure and the 3 
components of the velocity). For the magnetic field, however, 
we need to be more careful in order to enforce the divergence 
free condition. We use the zero-gradient condition only for the 
transverse magnetic field components B ± (perpendicular to the 
unit vector normal to the surface), and we use a linear extrapo- 
lation for the normal component B\\ (parallel to the unit vector). 
We make sure that the divergence of the magnetic field takes 
the same value on each side of the domain boundary, namely 
zero. 

This particular set of boundary conditions imposed on 
the magnetic field do not rigourously ensure that no electro- 
magnetic energy flux crosses the boundaries inwards (non- 
vanishing Poynting flux). We have checked that our results are 
not polluted by any artificial increase of the magnetic energy 
in the simulation coming from this choice of boundary con- 
ditions. The total magnetic energy increase during the whole 
simulation course coming from the boundaries is evaluated to 
be 10" 6 smaller than the magnetic energy increase measured in 
the galaxy at the end of the simulation. 

4. Disk formation with magnetic field 

When radiative cooling is turned on, our hydrostatic halo cools 
down from the inside out. Each shell looses pressure support 
and sinks towards the centre of the halo. The magnetic field 
lines are frozen in the free-falling plasma. Although initially 
the field lines are mainly vertical (perpendicular to the future 
galaxy), they are entrained and bent by the collapsing gas, so 
that the radial component finally dominates. Figure Q] compares 
the topology of the initial magnetic field and the magnetic field 
after 3 Gyr of evolution. One clearly sees that the field lines 
close to the disc are almost horizontal (with a strong radial 
component). The radial component changes sign above and be- 
low the plane, a natural consequence of the pinching of the field 
lines by the collapsing gas. Because of angular momentum con- 
servation, each collapsing shell will speed up, flatten and form a 
centrifugally supported disc. As we will show below, the galac- 
tic velocity field will quickly generate a toroidal field compo- 
nent. This toroidal field is a consequence of the field lines being 
wrapped up and amplified by the rotating disc. 

One should also note that the halo magnetic field keep its 
initial magnetic dipole-like structure, even after 3 Gyr evo- 
lution. This a direct consequence of the choice of our initial 
conditions. We will see that when a galactic wind develops, it 
strongly suppress the memory of the initial configuration of the 
field. 

We show in Figure [2] the initial magnetic field amplitude 
along the z axis of the halo. We consider in this study 3 different 
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Fig. 2. Magnetic amplitude along the z axis for different nor- 
malisations of the initial profile: #igm = 10" 6 yuG (dot-dashed), 
B lGM = 10 _5 yuG (dashed) and Bi GM = 10 -4 //G (dotted). The 
solid line shows for comparison the magnetic field that would 
have been in equipartition throughout the halo. 

cases: #igm = 10" 6 //G, #igm = 10" 5 yuG and #igm = 10~ 4 yuG. 
We also show in Figure [2] for comparison the magnetic field 
corresponding to equipartition in the initial halo. We define 
here the standard plasma parameter 



which measure the ratio of the thermal to the magnetic pres- 
sure. Equipartition means here J3 = 1. We see from figure [2] 
that for our 3 scenarios, the initial magnetic field is well below 
equipartition. 

After the disc forms, the magnetic field is amplified by 
rotation. If the resulting field is strong enough, the Lorentz 
force can significantly affect the dynamics of the flow, and al- 
ter the formation of the disc itself. The simulation with #igm = 
10" 4 yuG is quite extreme in this respect: the magnetic field is 
so strong that it prevents the formation of the thin, centrifugally 
supported disc, leaving instead a low density, highly magne- 
tised torus, in which star formation is very weak. The magnetic 
field lines become harder to deform by gas motions, except 
along the z axis, where field lines remain vertical and do not 
prevent gas from collapsing. This is clearly shown in Figure 
where we can compare the disc morphology in our 3 scenar- 
ios. The 3 maps of the J3 parameter are usefull to discriminate 
between the high magnetic field case #igm = 10~ 4 y^G, for 
which ft is very small in most of the torus region, demonstrat- 
ing that the torus is magnetically supported against gravity, and 
the low magnetic field case #igm = 10" 6 >uG, for which the 
disc we obtain is indistinguishable from the pure hydro case, 
and the magnetic field inside the disc is one order of magnitude 
below the equipartition value. Only in the intermediate case 
#igm = 10" 5 jdG do we get a thermally supported disc, with 
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a magnetic field roughly in equipartition only inside the disc 
08* 1). 

It is quite interesting to see thanks to these simulations that 
dwarf galaxies formation can be almost suppressed (or at least 
star formation within them) by a strong enough magnetic field 
in the IGM. It is even more interesting to notice that the value 
of the magnetic field that we require from cosmological simula- 
tion of galaxy cluster (#igm - 10" 5 - 10" 4 yuG) is precisely the 
same critical value above which dwarf galaxy formation can be 
regulated, if not suppressed, by magnetic fields. This point has 
important consequences that we will address in the conclusion. 

5. A Galactic Dynamo ? 

We now discuss more quantitatively the amplification mech- 
anism of magnetic energy in the galactic disc. W e use for that 
purpose the so-called ' 'Galac t ic Dynamo" theo ry dParkerll97l[ 



Wielebinski & Krausd fl993l: IShukurovl l2004l) . for which the 
mean large scale magnetic field evolution is governed by a 
modified induction equation 



dB 

- = Vx(hx5) + Vx (aB) + V x (77V x B) , 
ot 



(12) 



where parameter a represents the field amplification due to 
small scale turbulent motions driven by various small scale 
phenomena (supernova explosion, cl oud-cloud collision or col- 
lapsin g v ortex modes as explained i n lFerrierell992allEfstathiou 
200(JandEuIsrud & Zw eibel 2008) and parameter 77 represents 
magnetic diffusivity induced by turbulent transport and recon- 
nection. The standard a-£l dynamo theory relies on the a term 
to amplify the radial magnetic field B r and on the large-scale 
induction term to amplify the toroidal field B e . In the present 
simulation, however, we do not include any explicit a and 77 
terms in our induction equation. The original goal was to rely 
on our implementation of supernova feedback to induce heli- 
cal motions self-consistently. Probably because of our limited 
resolution (around 150 pc) and the fact that we use a subgrid 
model based on an effective EoS, we didn't get any a effect: 
field a mplification in our case is due only to differential ro- 



tation. iGressel et al.l (120081) . simulating a local patch a galaxy 
with much more resolution elements, obtained a galactic dy- 
namo driven by supernovae explosion and differential rotation. 
They outlined that the rotation is the critical driver for the dy- 
namo to operate an exponential growth of the field. 

An important issue is that we attempt to model a quiescent 
star-forming galaxy with an EoS mimicking the multiphase na- 
ture of the ISM gas, i.e. the gas pressure is artificially boosted 
to take into account the effect of supernovae heating of the 
ISM. The direct consequence is that the gas distribution is very 
smooth in the disc (no small scale clumping can develop), and 
as a consequence, supernova explosions have a very low im- 
pact on the ISM. A next step would be to simulate the evolu- 
tion of starburst galaxy, within which large clumps (~ 100 pc 
size) would be resolved and within which we could hope that 
supernovae are able to produce a small-scale dynamo effect. In 
order to get a global dynamo effect, additional physics would 
also be needed: non-ideal effects like Ohmic or turbulent dis- 
sipation and cosmic rays have been already identified as key 




Fig. 4. Total magnetic energy amplification as a function of 
time (a) for different values of the initial magnetic field (up- 
per panel) #igm = 10" 7 yuG (solid), #igm = 10" 6 yuG (dot- 
ted), Bigm = 10 _5 yuG (dashed) and Bi GM = 10 _4 /iG (dash- 
dotted), and (b) for different values of the Lagrangian refine- 
ment criterion with #igm = 10" 7 yuG (bottom panel): m res = mo 
is the reference run (solid line), m res = 5 x mo (dotted line), 
m res = 25 x mo (dashed line) and m res = 100 x mo (dash-dotted 
line). In each plot, the simple Q amplification model is over- 
plotted in red for Qq = 19 Gyr -1 as average angular velocity. 



ingredients to get a strong dynamo (see Hanasz et al ] l2009al for 
instance). 

An investigation of the magnetic energy growth within the 
disc reveals that there is no substantial increase of the magnetic 
field when SN explosions are allowed in the galaxy (figure 0. 
There is 3 effects: first, due to supernovae heating (both tur- 
bulent and thermal), the early magnetic field amplification ap- 
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Fig. 3. Slice through the (Oxz) plane of the gas density in log H/cm -3 (first row), of the temperature in log K (second row), of the 
magnetic amplitude in log yuG (third row) and the J3 parameter in log ft (fourth row) for an initial magnetic field of #igm = 10~ 4 //G 
(left column), #igm = 10" 5 juG (middle column) and #igm = 10" 6 yuG (right column) at t ^ 3 Gyr. Picture size is 40 kpc. 
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Fig. 5. Evolution of the magnetic energy (solid), thermal energy 
(dotted), and kinetic energy (dashed) for the simulation without 
star formation (black) and with SN feedback (red). The initial 
magnetic field is in both case #igm = 10~ 5 yuG, and energies are 
computed in a small cube of 15 kpc size around the galaxy. The 
blue solid line is the magnetic energy for the magnetic seeding 
simulation (see Appendix). 



pears weaker than in the no-SN case. As a result, the magnetic 
energy is slightly lower than without SNII. Second, the longer 
term magnetic field amplification appears slightly more effi- 
cient in presence of SN. Finally, the large-scale galactic wind 
carry some of the magnetic energy available in the disc, produc- 
ing a slight decrease in comparison of the simulation without 
feedback. As a result, the final magnetic energy value is almost 
equal to the final energy without SN explosions. This is a very 
weak and subtle effect. We notice that we reproduce roughly 
the same energy evolutio ns (kinetic, thermal and magnetic) 
than I Wang & Abell d2009l) . Interestingly, both kinetic and inter- 
nal energies are lower when feedback proceeds: this is partly 
due to the gas consumption by the star forming process, and to 
the ejection of some material outside from the galaxy. 

The evolution of the magnetic field is rapidly dominated by 
its tangential component defined in the cylindrical coordinates 
(r, and z, where z is the spin axis of the galactic disc). Thus 



d t B r ^0 and d t B z ^0 (13) 
when compared to the evolution of B e that rigorously writes 



d t Bg = - Bgd r v r - Bq8 z v z + B r d r VQ + B z 8 z vq 

- v r d r B e - —d e B e - v z d z B e -—B r . 
r r 



(14) 



For an axisymmetric field the amplification is done by the 
shearing terms, i.e. combining the third and last terms of this 
equation it leads to 



where we introduce the angular velocity Q(r). Then the radial 
shear governs the evolution of the magnetic field. We have ver- 
ified that the vertical shear component (rB z d z D) is negligible 
here. If the right-hand side terms of equation (TT3T) are not time- 
dependent, this equation resumes to 



B e = rB r d r £l x t . 



(16) 



If we assume that initially the magnetic field is dominated by 
the radial component, we can derive a very simple expression 
for the evolution of the total magnetic energy in the disc by 
integrating (B 2 + B 2 )l%n 



E B (t)/E B (0)=l + (£l G t) 2 



(17) 



d t Bg = rB r d r Q , 



(15) 



where =< (rd r Q) 2 > is the average (B 2 weighted) angular 
velocity squared of the galaxy. This simple analytical formula 
compares very well with the actual magnetic energy amplifi- 
cation we have measured in our simulations, as can be seen in 
Figure [4] Then one can infer that most of the energy comes 
from winding up the radial field within the disk. The best fit to 
our pure D amplification model on the magnetic energy evolu- 
tion (equation (fTTt ) from our simulations gives Do -19 Gyr -1 , 
which corresponds, using a simple dimensional analysis, to a 
disc magnetic scale length ^ 1.9 kpc for a circular veloc- 
ity of V200 = 35 km/s (with Qg - ^200/^)- This value is 
consistent with the angular velocity profiles measured at dif- 
ferent times through the course of the simulation (see figure [6]) . 
Only in our high magnetic field case, #igm = 10" 4 //G, do we 
see a significant deviation from this simple model. This is due 
to magnetic braking: angular momentum is removed from the 
disc b y Alfven waves propagating a long the more rigid field 
lines faennebelle & Fromangll2008 ). The angular velocity of 
the disc is therefore slowly decreasing. For all the other cases, 
magnetic braking is negligible and the magnetic amplification 
is very similar to the pure £1 amplification model. 

We have estimated the effect of mass resolution on the mag- 
netic amplification. Figured shows the measured amplification 
in 4 different cases: our fiducial with m res = mo, and 3 lower 
mass resolution with m res = 5, 25 and 100 x mo. The collaps- 
ing halo is therefore described with fewer and fewer AMR cells 
(3.10 6 , 8.10 5 , 4.10 5 and 3.10 5 total cells in decreasing order of 
m res ). The disc is however always refined at the maximum level 
of refinement, with a spatial resolution of 150 pc. As a conse- 
quence, only our lowest resolution simulation shows a signifi- 
cant deviation from the fiducial case, which can be considered 
as fully converged, within the framework of our sub-grid effec- 
tive EoS. 

A more detailed inspection of the magnetic field amplifica- 
tion shows that the contribution of pure compressional amplifi- 
cation of the magnetic energy (p 4/3 ) is dominant in the first 100 
Myr of galaxy formation due to the low velocity amplitude of 
the gas in the centre of the NFW profile. But after 200 Myr, 
this contribution to the magnetic field evolution, is negligible in 
comparison to the measured magnetic energy in the disc. So it 
is reasonable to not consider the amplification by compression 
over the whole galaxy evolution (which runs over Gyrs). 

Our simu lation results c ompar e very well with the work 
presented in K otarba et al. where these authors re- 

port a similar magnetic amplification using 2 different Smooth 
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phology is shown in Figure [71 We also represent in the same 
figure the magnetic pitch angle defined as 



Fig. 6. Mass-weighted average angular velocity in the galactic 
disc as a function of radius at t = 1 Gyr (solid line), t = 2 
Gyr (dotted line) and t = 3 Gyr (dashed line) for the simulation 
without star formation. 



Particle Hydrody amies (SPH) codes. They used various force 
softening length, down to 100 pc, quite similar to our grid spac- 
ing. They also used a similar physical model, since they didn't 
allow the gas to cool below 10 4 K. Although their SPH simu- 
lations were all initialized as an initial Milky-Way disc config- 
uration, their results are quite consistent with a pure £1 amplifi- 
cation. In a few cases only did they obtain a faster, exponential 
growth of magnetic energy. These cases all corresponds to nu- 
merical schemes where the magnetic divergence was allowed 
to deviate significantly from zero. They are probably associ- 
ated with numerical instabilities. In the majority of the SPH 
results, for which the magnetic divergence was well-behaved, 
the magnetic energy was growing like t 2 . 

Wang & Abel d2009h . on the other hand, have used an AMR 



code with a much better spatial resolution (down to 25 pc) and 
have included cooling processes down to 300 K. They were 
able to resolve, in a dwarf galaxy smilar to the one studied 
here, the fragmentation of gas clumps and the associated vor- 
tex modes. They however didn't consider star formation and 
supernova feedback. They concluded that the main amplifica- 
tion mechanism in their simulations was differential rotation: 
even at the 25 pc resolution they have reached, they have not 
captured any turbulent dynamo and their associated a effect: 
magnetic energy was also growing like t 2 in their case. 

To support even further the simple D amplification model to 
explain our simulation results, we now analyze in more details 
the magnetic field morphology in our disc galaxy. We consider 
for that purpose 3 different cases: 1- formation of the disc with 
only cooling, 2- formation of the disc with cooling and star for- 
mation and finally 3- formation of the disc with cooling, star 
formation and supernova feedback. The magnetic field mor- 



tan 0„ 



Be 



(18) 



P gives the angle between the magnetic field within the hori- 
zontal plane and its azimuthal component. P = 0° corresponds 
to pure azimuthal field lines and P = 90° to pure radial field 
lines. We must underline that pitch angles inferred from ob- 
servations are obtained with synthetic polarisation maps from 
radio emission. In the following, we extract a slice of the mag- 
netic field through the galactic plane to get the pitch angle as 
defined by equation ([T8K assuming that it provides a good com- 
parison with observational synthetic polarisation if the pitch is 
independent of the vertical height (it is roughly the case in these 
simulations). 

Without any star formation process, the magnetic field in 
the galactic disc is almost purely azimuthal: we can see from 
the left plot in Figure [7] that in this case the pitch angle is 
reaching a maximum value of P = -2°. The gas disc is sta- 
bilized by the effective EoS, so that no density perturbation can 
break the almost perfect azimuthal symmetry (only modes with 
m = emerge). After collapse, the radial component of the 
magnetic field points away from the centre of symmetry above 
the galactic plane and towards the centre below the plane. 
Using Equation (fT5K we see that the toroidal magnetic field 
also changes sign above and below the plane: the field is said to 
be anti- symmetric. This kind of particular symmetry is called 
the AO mode (anti- symmetric and m = 0) and emerges natu- 
rally in halo without strong dis c-halo integrations (like galac - 
tic winds or galactic fountains) (IBeck et al.lll996clBeckll2009l) . 
The Milky- Way, for example, is believe to be a good exam 



pie of such a A O magnetic field in the halo (IHan et al 
Sun et al.lEQ08h . 



1997 



The situation appears to be quite different when star for- 
mation in taken into account. Because mass is removed from 
the gas disc to form stars, our effective EoS reaches lower 
density and the gas sound speed is smaller. The disc Toomre 
parameter get smaller, closer to 1, and a strong spiral mode 
can develop. We see from the central plot of Figure [7] that the 
magnetic field is not axisymmetric anymore. The strong den- 
sity perturbation triggers fluctuation in the angular velocity, 
so that D(r) is not monotone anymore, but fluctuates around 
V20o/ r We can see again from Equation (TT3T) that, when B r 
changes sign, the corresponding toroidal component can de- 
crease, or even change sign. We see indeed in Figure |7] that 
the toroidal field almost reverts its main direction between the 
spiral arm and the inter- arm regions. As a consequence, the 
pitch angle amplitude \0 P \ is rather high in the low density, inter- 
arm regions \0 P \ ^ 40° (with both negative and positive signs) 
and is lower within the spiral arms with \0 P \ - 20° (always 
with negative values). These values are in good agreement 
with the observed pitch angles in nearby galax ies (IBeck et al 



with t he observed pitch i 
1996ilRohdeeTaT]ll999[ 



Patrikeev et al. I l2006h . The same re- 



lation between the pitch angle and the spiral arms is also ob- 
served within spiral galaxies: magnetic field is aligned with the 



optical spiral arm s (Krauseetal. 19891: iNeininger et al.lll991 



Berkhuijsen et al ]l997t lBeckll2007l) . and it also appears some- 
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Fig. 7. Magnetic amplitude in log jiG with magnetic vectors overplotted (upper panels) and the pitch angle 6 P in degrees (bottom 
panels) for the galaxy without star formation at t ^ 3 Gyr (left panels), for the galaxy with star formation and without feedback 
at t = 4 Gyr (middle panels), and for the galaxy with star formation and supernova feedback at t = 3 Gyr (right panels) in the 
upper part of the equatorial plan. Initial magnetic field is #igm = 10" 5 yuG. Picture size is 40 kpc. 



times that stronger ordered magnet i c fields are seen w ithin the 
inter-arm regions dEhle et al.lll996l:lFrick et al.ll2000h . 



The figure |7] clearly shows positive and negative signs of 
the pitch angle. Therefore a simple spiral density flow can not 
explain the observations, which show only negative pitch an- 
gles. But spirals may contribute to the increase of the pitch an- 
gle compared to a pure axisymmetric rotation, at least at the 
inner edge of the spirals. Thus, the always negative pitch angle 
seen in the observations is one of the arguments favouring the 
dynamo. 

The Magneto Rotational Instability (MRI) could also 
be a growth mode of the magnetic field in galaxies 



(IKitchatinov & Riidiger 2004) or to sustain turbulence in the 
disc ( Piontek & Ostrikei 120051) . The MRI can give a coherent 
picture for the pitch angle of polarised emission if the mag- 
netic field has sev eral reversals along the vertical direction 
(Elstn er et al ] l2009h . The MRI needs a strong toroidal field to 
be efficient, and it is clearly shown that the field in our galaxy 
simulations is preferentially poloidal at late times (see figure[7]). 
However there is a possible MRI channel for the magnetic field 
growth in the early evolution of the galaxy (< 300 Myr) when 
the circular magnetic component is still lower or comparable to 
the vertical field, then, later on, the amplification by differential 
rotation proceeds. But it is beyond the scope of this paper to 
disentangle this early MRI from the pure differential rotation 



suggested before and would require simulations with higher 
resolution to properly describe the MRI driven- turbulence. 

When supernova explosions are finally included in the 
model, the previous magnetic field topology is conserved. The 
spiral arm is somewhat weaker, because the increased veloc- 
ity dispersion due to supernova blast waves stabilizes the disc. 
The disc appears less coherent and more perturbed by small 
scales perturbations. Note that this is precisely the small scale 
turbulence we need for the a effect, but despite this small 
scale perturbations, the radial component did not grow sig- 
nificantly in our simulation. The magnetic energy evolution 
conserves the same properties than in the run without star 
formation: the differential rotation of the disc is the main 
driver of the galactic magnetic field in these simulations. We 
probably need sub-parsec scale resolution in order to resolve 
supernova induced c yclonic motions that can drive a turbu- 
lent dynamo, like in iBalsara et al.1 ( 2004 ). Besides the lack 
of small scale amplification of the radial component of the 
field, our simulation recovers most of the qualitative features 
of Galactic dynamo models, that predict similar pitch angle dis 



tribut i on and magnetic field topologies (IKr asheninnikov a et al. 
19891: bonner & Brandenburg! 1 1990: Elst ner et al.lll992h . The 



existence of this small scale a effect in galaxies was re- 
cently shown by direct simulations of iGressel et al.1 (12008b 
but the problem of its saturation is still under debate 



12 



Dubois and Teyssier: Magnetised winds in dwarf galaxies 



53 



51 



50 



I 


■ i i i 


: ! 





t (Gyr) 



Fig. 9. Magnetic energy flux as a function of time for #igm = 
10" 5 //G measured at a radius of 4r s (dotted line) and 9r s 
(d ashed line) from the halo centre. The analytical prediction 
of iBertone et al.1 d2006h is overplotted as the red solid line. We 
have overploted the result for the direct magnetic energy injec- 
tion with zero initial magnetic field (solid line) as described in 
Appendix lAl 



( Brandenburg & Subramanianl 120051 for example). The value 
of the parameter a needs to be rather high in order to obtain 
a fast enough galactic dynamo. The conservation of the total 
magnetic helicity in the galaxy makes large value for a rather 
difficult to obtain if there is no exchange of magnetic helic- 
ity with the halo. In this paper, we rather focus on the other 
consequence of supernova feedback, namely the formation of 
a strong galactic wind, and how this wind transports magnetic 
energy. 

6. Magnetic field in the galactic wind 

We now study in more details the properties of the supernova- 
driven galactic wind that forms during our simulations. We re- 
strict ourselves to the analysis of the case #igm = 10" 5 //G: 
we obtained in this case a normal star-forming disc, and in the 
same time a magnetic field close to equipartition in the gaseous 
disc. A larger value of #igm would result in the formation of 
a non-star-forming torus, while a smaller value would give 
rise to a weakly magnetised galactic disc, significantly below 
equipartition. 

In agreement with the simulations shown in 



Dubois & Tevssierl d2008bl) . the galactic wind appears at 
1.5 Gyr and fully develops during the next Gyr. Figure [8] 
shows 4 different snapshots (from 1 to 3 Gyr) of the magnetic 
field in a plane perpendicular to the disc. These snapshots 
illustrate the 3 di fferent phases of the galactic wind, described 
in more details in Dubois & Tevssierl d2008bl) : 1- the bubbling 
phase, when the supernova luminosity is too weak to break the 
ram-pressure of the inf ailing gas it = 1 Gyr in Fig. [8]), 2- the 



snow-plow phase, when the supernova luminosity is greater 
than the ram-pressure of the infalling gas, and the wind starts 
to propagate through the halo (t = 1.5 Gyr and t = 2 Gyr in 
Fig. [8]) and 3- the blow away phase, when the wind has reached 
its final, stationnary shape with a characteristic noozle-like 
structure (t = 3 Gyr in Fig. [8]). 

We can see from Figure [8] that the magnetic field in the 
wind is highly turbulent, with an average amplitude around 
10" 3 yuG, but showing strong variations within the wind, with 
highly magnetised filament (as high as 10 _1 yuG) entrained by 
low density bubbles (as low as 10" 4 yuG). The turbulent struc- 
ture of the magnetic field is typical of convective flows, with 
buyoantly rising vortices and their associated shearing motions. 
The magnetic field lines appear however mainly colinear to 
the flow, a magnetic co nfiguration that is actually observed in 
galax i es with outflows (Brandenburg et alJ 1 1993 ; IChyzy et aT 



20061: iHeesen et al.ll2009h . Whether the turbulent flow we ob 



serve in the galactic wind amplifies the magnetic field is un- 
clear in our simulation. As we discuss below, the wind mag- 
netic energy is clearly decaying because of the overall expand- 
ing flow. Moreover, in the wind region, our spatial resolution 
is slightly degraded (around 300 pc), because the grid has been 
derefined in this low density region. This suggest that we are 
not resolving the small scale magnetic energy amplification due 
to the turbulent convective flow, leading to a potential underes- 
timation of the magnetic energy in the wind. 

In order to estimate the amount of magnetic energy that the 
wind can extract from the magnetised galactic disc, we have 
measured the magnetic energy flux in 2 thin shells of radius 
4r s 20 kpc and 9r s 45 kpc. Results are shown in Figure [9) 
After the wind starts (around t = 1.5 Gyr), the magnetic en- 
ergy flux quickly jumps to a value of 10 52 erg/Gyr, and stays 
roughly constant during the next Gyr. The magnetic energy 
flux is decreasing when going to larger radii (roughly as r" 1 ), 
demonstrating that the effect of expansion is overcoming field 
amplification in the wind turbulence. It seems that the flux still 
grows after 2 Gyr. It is possible that when the wind reaches a 
stationary regime, this growth is correlated with the increase of 
magnetic energy within the disc (see figure [5]). 

It is enlightening to interpret our n umerical result s in th e 
framework of the analytical theory of IBertone et al.l d2006h . 
This will also guide our discus sion on magnetic fie ld amplifi- 
cation in general. According to IBertone et al. 1 d2006h . the mag- 
netic energy flux which is injected at the base of the galactic 
wind can be written as 



-B,in 



B 

on 



(19) 



where B m is the injected magnetic field, Rq - 7 kpc is the 
galactic di sc radius and is th e wind velocity. In their analyti- 
cal model, Bertone et al. I d2006h assumed first that the magnetic 
energy could be considered as some isotropic pressure term. 
We see from the previous discussion that the convective turbu- 
lence in the wind justifies this assumption. They assume also 
that the injected magnetic field scales with the injected wind 
gas density as 



B\n — Be 



PinV 
PG) 
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(20) 
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Fig. 8. Magnetic field in log //G units for the galaxy with star formation and supernova feedback with #igm = 10" 5 y^G at different 
epochs t = 1 Gyr (upper left), t = 1.5 Gyr (upper right), £ = 2 Gyr (bottom left) and t = 3 Gyr (bottom right). Picture size is 40 
kpc. 



where Bq and po are respectively the typical magnetic field 
and gas density inside the galactic disc. Finally, they related 
the wind gas density to the wind mass flux by 



M w = p in 47rR G u w . 



(21) 



Using the 3 previous equations, we can predict the magnetic 
energy flux in the wind to be 



Ebm = £b,< 



4u w I M W R G 



Hq \3u w M] 



ISM 



(22) 



where Hq ^ 300 pc is the disc thickness and Eb,g is the disc 
magnetic energy. We have plotted the predicted magnetic en- 
ergy flux corresponding to Equation (l22t in Figure [£l using the 
measured mass flux (M w 0.035 M©/yr), disc magnetic energy 



(E b ,g - 5 x 10 53 erg) and disc mass (Mi S m - 10 9 M ). We see 
that this formula is quite accurate, given the number of sim- 
plifying assumptions we have made. It slightly underestimates 
the measured flux, by less that a factor of 2. If we integrate this 
magnetic energy flux over ^ w - 1 Gyr of wind activity, we find 
that roughly 6 % of the disc magnetic energy can be extracted 
by the galactic wind and funneled into the IGM. We note 6 W this 



global wind efficiency. In the analytical model of lBertone et al 
d2006l) . this efficiency writes 



6w = 



Hq \ 3w w Mism / 



(23) 



A wind mass injection of 0.03 M©/yr is a rather low value com- 
pared to starburst galaxies where wind are usually observed. 



14 



Dubois and Teyssier: Magnetised winds in dwarf galaxies 








1 I 1 1 : 


-1 


[ IGM 


HALO i 














« -3 

o 






-4 


||J 




-5 




1 . i .... i ... i ... : 



-2.0 



-2.5 



-3.0 



-3.5 



-4.5 



log p/(p) 




2 4 
log p/{p) 



-2.51 



-3.0 



-3.5 



-4.0 



-4.5 



-5.01 



Fig. 10. Magnetic amplitude versus gas over-density diagram at 
t = Gyr (upper panel) and t = 3 Gyr (bottom pannel) for the 
galaxy with star formation, supernova feedback and #igm = 
10" 5 yuG. Black solid line is for the initial p 2/3 law and black 
dashed line is for the initial p 2/3 law amplified by a factor 60. 



Our simulation can be considered as a quiescent wind case. 
Typical starburst- driven winds show strong mass outflow as 
large as 10 M Q /yr dMartinll 1 99 8l 1 1 999h . On the other hand, the 
burst duration is pro bably much smaller, with t w 100 Myr 
(Ide Grijs et al ] l200lh . Overall, the efficiency of magnetic en- 
ergy extraction probably lies above e w ^ 10 % for a typical 
star forming dwarf galaxy. The other extreme case one could 
also consider is when the starburst is so violent that the whole 
galaxy is destroyed. This scenario is likely to happen in the 
early universe inside small primordi al galaxies, as shown by th e 
recent cosmological simulations of iMashchenko et aD d2008 ). 
In this case, the efficiency obviously reaches 6 W = 100 %. 

We now want to estimate the magnetic field in the IGM 
that will result from the enrichment of the galactic wind. 
For this purpose, w e use the analytical model developed by 
Bertone et al. I d2006h and used recently in a cosmological sim- 
ulation by iDonnert et al. I d2009h . We write for the evolution of 



the magnetic energy in the wind blown bubble the following 
equation 



dt 



Z7 ^ W Z7 

^B,in - — : 



(24) 



£b,ui is the magnetic energy injection rate, measured at the base 
of the expanding wind and R w is the bubble radius (R w its time 
derivative). The second term on the right-hand side of the pre- 
vious equation stands for magnetic energy losses due to the 
expansion of the bubble. This equation also assumes that we 
can neglect the turbulent amplification of the magnetic field, 
which turned out to b e true in the s i mulati ons we have pe- 
formed here. Following IDonnert et alJ d2009h . we consider the 
rather pessimistic case of a fast bubble expansion and writes 
R w /R w = I It. We also consider the solution of Equation (l24l) 
for which the initial bubble energy vanishes 



^B,in - 



-£b,g • 



(25) 



We note that, in this scenario, half of the injected energy has 
been lost due to the expan sion of the field lines . Using the cos- 
mological wind model of iBertone et al. I d2006h . we consider a 
final wind radius of R w = 800 kpc and a wind filling factor of 
f w = 100 % (both typical values for the universe at redshift 
zero), so that the final magnetic field in the IGM writes 



l3e w E] 



B,G 



Rl 



1.1 x htVg. 



(26) 



In conclusion, the galactic wind has enriched the IGM with a 
magnetic field more than one order of magnitude larger than its 
initial value. Interestingly enough, the wind bubble magnetic 
field is proportional to the disc magnetic field. Indeed, the disc 
magnetic energy depends on the toroidal magnetic field in the 
disc as follows 



Eb,( 



8tt 



ttR 2 g H g . 



so that the new IGM value now writes 



B 



IGM 



\?)€ w R 2 c Hq 



(27) 



(28) 



This last formula stands as a nice and concise summary of our 
work: we have shown with detailed MHD simulations that the 
final magnetic field in the expanding wind bubble inherits its 
value from the parent dwarf galaxy, where an £1 effect has am- 
plified the field by almost 2 orders of magnitude from its initial 
value. 

7. Discussion 

Using a rather idealized set up to model dwarf galaxy forma- 
tion with MHD, we have followed the evolution of the mag- 
netic field, from the collapse of an equilibrium halo, to the for- 
mation of a galactic disc and finally to the ejection of magnetic 
field lines by a supernova driven galactic wind. In order to illus- 
trate further the mechanism powering up this cosmic cycle, for 
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which dwarf galaxies play the central role, we have plotted in 
Figure [10| the magnetic phase space diagram, showing a mass 
weighted histogram as a function of magnetic field amplitude 
and gas overdensity. First, when the halo forms, the magnetic 
field is amplified by gravitational contraction from its initial 
IGM value to its final disc value 



B r = Biqm 



/ \2/3 

— 

\PlGM / 



(29) 



We use here only the radial component B r , because just after 
collapse, the magnetic field is essentially radial in the disc (see 
discussion above and Fig. []]). We can see in the upper part of 
Figure [10| our initial configuration, where B and p are follow- 
ing this tight relation. Note here that we have neglected the 
additional amplification of the magnetic field due to turbulent 
motions in the halo. If we note R the initial Lagrangian radius 
of the collapsing halo and A its spin parameter, we k now from 



the sta ndard cosmological disc formation model of iMo et al 
dl998l) that the final disc radius writes 



200 1 / 3 



so that the initial disc radial component writes 

200 2/3 



B r ^ B\Q M - 



A 2 



(30) 



(31) 



In the lower panel of Figure [TO1 we see that the toroidal field 
in the disc has grown linearly with time due to differential ro- 
tation (the D amplification), by a factor of 60, corresponding to 
the growth rate (Eq. [16]) integrated over 3 Gyr. We finally get 
the growth rate of the toroidal magnetic field in the disc as a 
function of the initial magnetic field in the IGM as 



200 2/3 

d t B e - B r £ls ^ Biqm — — • 



(32) 



In the lower panel of Figure[lOl we see that this amplified mag- 
netic field is funneled out of the galactic disc, following a p 2/3 
relation back to the IGM. This is the "galactic wind" branch 
that ends up the overall cycle. We have seen that the IGM mag- 
netic field has been amplified by one order of magnitude during 
this cycle of roughly 3 Gyr. 

Although our simulations, in conj unction with the wor k 
presented in iBertone et al. I J2006I) and Ponnert et al. d2009l) . 
provide a nice explanation for the origin of an IGM magnetic 
field, at a level as high as 10" 4 yuG, assuming equipartition mag- 
netic fields in dwarf galaxies, we still need to explain the origin 
of these intense galactic fields, given that we start our cosmic 
evolution wih a very small initial value B\gm - 10~ 14 yuG just 
after reionization. We can't rely on the £1 amplification alone, 
since in this case, the magnetic field grows only linearly with 
time. In the context of isolated galaxies, the only source of ad- 
ditional amplification comes from small scales turbulent and 
cyclonic motions (the a effect). This turbulent dynamo is be- 
lieved to amplify the radial component of the galactic field, the 
toroidal component being generated and amplified by the Q ef- 
fect. 

We know however that galaxies are far from being iso- 
lated. They form by accretion of diffuse IGM gas, usually 



in the form of den se 
dDekel et alJ l2009h 



cold filaments known as cold streams 
but also by accretion of satellite galax- 
ies, orbiting around the central galaxy before being incorpo- 
rated more or less violently into the final disc. These 2 impor- 
tant pro cesses define the hierarchical scenario of galaxy for- 
mation dSilk & White! Il978l: Iwhite & Reesl Il978l) . and prob- 
ably explain most of the properties of p resent day galaxies 
dGovernato et al.ll2007l: iMaver et alJl2008h . 

We propose here to look at the galactic dynamo problem 
from a different angle. Instead on relying only on the small 
scale dynamo to amplify the radial component of the field, we 
should also consider as a source of radial field: 1-the toroidal 
magnetic field coming from accreted satellite galaxies and 2- 
the accreted IGM magnetic field previously enriched by galac- 
tic winds or 3- the accreted IGM magnetic field previously en- 
riched by entirely destroyed progenitor galaxies. On one hand, 
we know from Equation (f28l) that the enriched magnetic field 
is proportional to the toroidal field of the parent galaxy, and on 
the other hand, we know from Equation d32b that the growth 
rate of the toroidal field of the main galaxy is proportional to 
the initial IGM magnetic field. 

Although a correct description of satellites and cold streams 
accretion requires detailed cosmological simulations of galaxy 
formation with MHD, we can speculate here that our proposed 
cycle of accretion and ejection of magnetic field lines (see 
Fig. [TO]) probably provides a dynamo loop, for which the radial 
field component is amplified by accretion of satellite galaxies 
and of the magnetically enriched IGM. A very rough and spec- 
ulative analytical description of this effect can be obtained by 
injecting Equation (f32l) into Equation 1 



dtBiGM - Biqm - 



3e w R 2 G H G R 2 



*g 



^g-^igmV^w^g, (33) 



where we have exploited the result of IBertone et al.l d2005l) that 
R w ^ 3R (a reasonably good approximation according to their 
Figure 2) and we have assumed a typical disc aspect ratio of 
Hg/Rg - 1/20 and a typical spin parameter A ^ 0.05. We 
see from the previous equation that now the magnetic field am- 
plification is exponential, and not linear anymore, because we 
have connected the radial component to the progenitor galax- 
ies' toroidal component. Altough we have clearly not demon- 
strated that this cosmic dynamo actually works, we want to 
stress that the cosmological environment might play an impor- 
tant role, namely by accretion of amplified field lines from dif- 
fuse gas or satellite galaxies. The central engines of this cos- 
mic dynamo are the dwarf galaxies, because they are easily 
disrupted by galactic winds, and therefore provide a important 
source of diffuse radial field, and also because they are the pro- 
genitors of our Milky Way Galaxy. 

Some caveats to this Cosmic Dynamo scenario must be 
pointed out. We assumed that #igm is regenerated in a pure 
radial mode to get that exponential growth. Even though mag- 
netic field lines in large-scal e bipolar outflows a re mainly per- 
pendicular to the disc (see iHeesen et al.l 120091 for example), 
they carry turbulence that produces non-radial component of 
the field. Propagation of the outflow can be affected by the par- 
ticular geometry of the galactic environment and leads to irreg- 
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ular distributions of the magnetic field components. Moreover 
magnetic fields can contribute constructively or destructively 
to the Biqm when large-scale bubbles percolate and can sub- 
stantially alter the final configuration of the field. A very rough 
way of taking into account these effects on the radial compo- 
nent is to consider that each component of the field is equally 
distributed, then equation (l33t is simply affected by a factor 
1/V3. 

While the a-£l dynamo comes from the connection be- 
tween small scale turbulence and galactic differential rotation, 
the proposed Cosmic Dynamo originates from the coupling of 
a large scale gas cycle (cosmological accretion and wind ejec- 
tion) and differential rotation. We see from Equation (f33l) that 
this effect is probably very fast, especially at high redshift, for 
which the angular velocity of dwarf galaxies is much higher 
than today (Qg °c (1 + z) 3/2 ). The actual growth rate will how- 
ever depend crucially on many complicated aspects of galaxy 
formation that are beyond the scope of this p aper, such as 



the w ind filling factor as a function of redshift (IBertone et al 



20051) . the field geometry during satellite mergers, the exact 
3 -dimensional propagation of wind bubbles within the cosmic 
web, etc. Here it is a priori assumed, that the galactic wind ex- 
pansion in conjunction with some mergers happen within the 
3 Gyr. Thus it has to be shown, that a similar growth could be 
achieved with an initial field configuration, which is similar to 
the evolved wind field. 

We believe that this Cosmic Dynamo, in conjunction with 
a strong galactic dynamo, should be able to account for the 
fast magnetic field amplification we see in our universe, espe- 
cially in light of the recent meas urements performe d on high 
redshift galaxies up to redshift 1 (IBernet et al. Because 
the existence of an efficient a-£l Dynamo is questionable in 
dwarf galaxies due t o low rotati on of the disc and the lack of 
relativistic particles ( Kleidfl991 ). it is important to explore al- 
ternative scena rios. It is possible to rely on the ideas discussed 
by lReesI ( 1987b , for which stars are able to generate their own 
magnetic field by stellar dynamos. We have tested this scenario 
in the Appendix, showing that for reasonable assumptions, we 
can reach magnetic field close to equipartition in a few Gyr. 

We have also demonstrated with our cooling halo simula- 
tions that if the IGM magnetic field is high enough (namely 
#igm - 10~ 4 //G), the Lorentz force can prevent the forma- 
tion of a new generation of dwarf galaxies. The formation of a 
magnetically supported torus prevent any subsequent star for- 
mation, and therefore breaks the central engine of our Cosmic 
Dynamo. This has the interesting consequence of providing a 
saturation mechanism, explaining naturally why the Cosmic 
Dynamo leads to a final value of #igm - 10~ 5 to 10" 4 yuG. 
While the a-£l dynamo saturates when the small scale turbu- 
lent magnetic field reaches equipartition in the disc and mag- 
netic teiision_be^mes^ 

(see ISubramanianl 1 19981 iBrandenburg & Subramanianl 2005 ). 
the Cosmic Dynamo saturates when the large scale, cosmolog- 
ical magnetic field prevents dwarf galaxies' disc formation, or 
severely affects star formation within them. 

The evolution of magnetic fields in galaxies can now be ad- 
dressed in a cosmological context. In a more or less near future, 
we propose to extend the present work using much higher reso- 



lution simulations, the ultimate goal being to resolve the small 
scale a effect in our isolated, cooling halo simulations. Another 
ambitious objective would be to address the same problem in 
a fully cosmological environment, capturing the a effect with 
some subgrid model, but accounting properly for satellite ac- 
cretion and wind bubble propagation. 
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Appendix A: Stellar magnetic feedback scenario 

iReesI dl987l) was the first to suggest that Biermann battery ef- 
fects at the surface of stars could create magnetic fields from a 
zero field. We also know that co nvective stars can amplify the 
magnetic energy extremely fast (Brun et al. I l2004 close to its 
equipartition value by dynamo effects like in the Sun. In this 
scenario, massive stars > 10M© which ones will release their 
magnetic energy in the ISM in powerful explosions, are the 
key ingredient for generating galactic magnetic fields. In the 
most massive stars, it seems that there is no convective enve- 
lope allowing for a very fas t dynamo, but differ ential rotation 
could be extremely strong (IMaeder et al.ll2008l) and can sub- 
stantially increase any seed field. In so me massive stars the con- 
vective zone can be in the c entral p art (Zah n et al.1 120071) and the 
dynamo can be effective (iCharbonneau & MacGregorll200ll) . 
These stars have also a very short life-time (a few Myr), any 
amplification process must have an ^-folding time extremely 
small. Some models predict that instabilities in the external ra- 
diative en velope are able to sustain a sufficiently fast growth of 
the fi eld dMacDonald & Mullanl Eool iMullan & MacDonald 
20051) . Though it remains difficult to conclude whether mag- 
netic field are appreciably amplified, there are increasing ob- 
servational evidence for strong magnetic fields in massive stars 
dAlecian et al.ll2008l:lHubrig et al.ll2008k . 

When a supernova is exploding, its magnetic field is diluted 
in a larger volume. We can approximate its final value by con- 
siderating the flux conservation of the magnetic field. The uni- 
form and averaged magnetic field in a spherical transformation 
must follow 



HI* 



(A.l) 



For example the Crab nebula with B\ ^ 100//G within a 
R\ ^3.4 pc radius (Ke nnel & Coronitil[l984]) will have a Bi = 
5.10" 2 yuG field within the radius R2 = 150 pc of a superbubble. 
The case of the Crab nebula is a very strong particular case: a 
fast rotating pulsar inside the centre of the remnant could be 
responsible for this very high magnetic field. It is particularly 
difficult to ex plain how strong ma gnetic fields B\ ~ 10 9 yuG in 



massive stars dAlecian et al.ll2008L R\ ~ 5R Q ~ 3.10 6 km) can 
produce so high values in supernova remants (B2 ~ 10~ 5 //G 
for R2 ~ 1 pc). This prediction does not take into account 
the possible amplifications stages during the rapid expansion of 



the blast wave, magnet ic field can grow due to the cosmic ray 
pressure ( Parker 1992|), some strong B iermann b attery effect 
dHanavama et alJl2005h . or turbulence dBellll2004 that would 
allow stronger values in supernova remnants. 

The basic idea to implement the release of magnetic field 
by superbubbles coming from supernova explosions is to put 
a fixed total magnetic energy within the bubble (ensuring that 
V.B = 0) with a random direction. We assume that within a 
^bubble - 150 pc radius, each supernova contributes construc- 
tively for a certain amount of magnetic energy to the final re- 
lease. We arbitrary choose that each superbubble (which cor- 
responds to an assembly of several individual supernova ex- 
plosions) releases a 10" 3 yuG field following a robust numer- 
ical implementation. This particular choice of 10" 3 yuG mag- 
netic seeding is entirely heuristic and is set to match jiG field 
after 1 Gyr of galactic evolution. To get a rough idea of the 
magnetic release by a single supernova inflating a larger bub- 
ble, we have to know the number of supernovae per star cluster 
particle Nsn,* = m* x ^/sn/^sn - 100, where we have as- 
sumed that M S n = 10M o and for ra* = n x Ax 3 ^ 1O 4 M . 
To simplify if we assume that each single supernova produces 
a magnetic field with the exact same intensity and direction, it 
comes that the associated magnetic field per supernova (10 M ) 
is #sn - 10" 5 yuG (in a 150 pc radius bubble). We do not take 
into account the destructive contributions (magnetic field in op- 
posite directions) but we do not also take into account the sub- 
stantial increase coming from shocks between supernova ex- 
plosions that generate instabilities, shear flows, cosmic rays, 
etc. 

To implement a method of direct injection of the magnetic 
energy in the code, we must find a form of the magnetic field 
that ensures V.B = 0. For that purpose we will design a dipolar 
moment configuration such that the resulting magnetic field is 
embedded in the radius of the supernova explosion, and that the 
magnetic energy is not divergent in the centre in order to pre- 
serve the stability of the scheme even at different resolutions. 
Let us recall that the magnetic potential vector created by a 
unique magnetic dipole follows 



A(r) = 



4n 

yu r 3 



-(m A r) , 



(A.2) 



where m is the magnetic dipolar moment. If we take m as a 
constant value, one can easily see that magnetic field will be 
infinite in the centre, then we must choose a certain distribu- 
tion of dipolar moment elements such that the overall resulting 
moment respects the compact and convergent criterions. For 
that purpose we need m to grow as fast as r 4 and that sharply 
decreases down to zero at a certain radius tsn 



m = r exp 



rsNj 



m . 



(A.3) 



With such a distribution of m, we can compute the magnetic 
field as 



B = V A A. 



(A.4) 



The resulting magnetic field computed on a 512 3 grid within 
a 5rsN box size (resolution element is 0.01 tsn) is plotted on 
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Figure IA.U We can renormalize the dipolar moment distribu- 
tion such that the magnetic field gives a unit magnetic energy 
in regards to the measured magnetic energy E m , here 



4n 



m 



f 6fi E ni 



0.85. 



(A.5) 



This value is not very sensitive to the numerical resolution, for 
example, if Ax = 0.3rsN, then mo = 1. That means that even if 
we poorly resolve the bubble (it is the case in our simulations 
tsn - 2-3 Ax), the final magnetic energy is very close to the 
converged one. 

We also add a saturation factor to avoid that the magnetic 
field excessively grows above the equipartition value, such that 
the resulting field is multiplied by 



fsat — 



1 



1 + 



(A.6) 



where B is the value of the magnetic field inside the bubble and 
B Qq is the magnetic amplitude corresponding to its equipartition 
state with internal energy. For each buble, the / sat factor is the 
maximum of all individual / sat ,i factors computed on each cell 
within the supernova radius. Such saturati on effects have been 
outlin ed in galactic dynamo simulations dCattaneo & Hughesl 
1996h . 

We produced the exact same simulation than in Section [6] 
with no initial magnetic field. Magnetic fields are randomly 
generated in the disc but they tend to align with the global ro- 
tation within the disc as one can see in Figure lAl2l After 3 
Gyr, the amplitude has reached a value of a few ~ 10~ 2 fiG 
in the central parts and its configuration is no more like a AO 
type. As there is no more initial magnetic field, the AO con- 
figuration does not appear. Moreover the randomisation of the 
magnetic field injection destroys any particular symmetry and 
reconnects the field through the equatorial plane. This effect 
explains the smaller scale of the magnetic field fluctuations but 
the global differential rotation of the galaxy maintains a pre- 
ferred component of the field along the rotating vector. A di- 
rect consequence of the differential rotation is the amplifica- 
tion applied on magnetic seeds. A computation of the injected 
magnetic energy from supernovae with the star formation rate, 
gives E{ n j ^ 10 48 erg /Gyr. Though the magnetic energy in the 
disc should reach 10 48 erg after 1 Gyr, figure \5\ shows that the 
energy measured in the galaxy for this simulation (blue line) 
is ~ 3.10 49 erg after 1 Gyr. This is an evidence that the dif- 
ferential rotation is still amplifying the seed field generated by 
supernovae. 

Radial contribution is non-negligeable, the pitch angle \6 P \ 
is everywhere greater than 10° and goes up to 90° in some re- 
gions of the galaxy. The 10° minimum value is predictable: 
if the mean field amplitude is about ~ 10" 2 yuG preferentially 
toroidal, putting a 10" 3 yuG field in the supernova remnant along 
the radial component produces a \6 P \ ^ 5° pitch angle. Then 
the successive explosions and the differential rotation rise this 
value up to 90°. 

As shown in figure [9j the magnetic energy flux is lower 
than when the simulation starts with an initial magnetic field 




Fig.A.l. Magnetic amplitude for the (IA.3b dipolar distribu- 
tion where m is aligned with the zaxis. Upper panel is a slice 
through the (Oyz) plane and bottom panel is through the (Oxy) 
plane. Magnetic vectors are overplotted. Picture size is 5rsN- 



by a factor ~ 40, but this value depends a lot on the partic- 
ular choice of the magnetic seeding B$n = 10" 5 //G. Thus, 
with a EB,m - 10 51 erg. Gyr -1 at the bottom of the wind, we 
can predict that the final magnetic field within the hot tenu- 
ous bubble will reach ~ 10" 5 yuG. It corresponds to the initial 
cosmol ogical seeds required in large scale cosmolog ical simu- 
lations (Dolag et al. 2005; Dubois & Tevssiedl2008al) . Relying 
on the assumption that some supernova remnants (Crab neb- 
ula, Vela nebula) are permeated with strong magnetic fields 
dKennel & Coronitilll984 iHelfand etaPl200lk we verify that 
this mechanism can magnetise a dwarf galaxy and therefore 
produce a magnetised wind to enrich the IGM. Using a very 
low magnetic seed within supernova bubbles (5000 weaker 
than the field within the Crab nebulae), we showed that the 
galaxy rapidly reaches 1 fiG in ~ 1 Gyr and is able to fill 
the IGM with a large-scale magnetic field of ~ 10" 5 yuG. It 
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is therefore sufficient to explain an early magnetisation of the 
Universe. 




Fig. A.2. Magnetic amplitude in log fiG units with magnetic 
vectors for the galaxy with star formation, supernova feedback, 
Bigm = and a magnetic seeding #sn = 10" 5 yuG at t ^ 3 Gyr 
in the (Oyz) plane (upper panel), in the (Oxy) plane (middle 
panel) and the pitch angle (bottom panel) . Picture size is 40 
kpc. 



